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Abstract 

In this work, based on quantum operator Hermite polynomials and Weyl’s mapping rule, we find a 
generation function of the two-variable Hermite polynomials. And then, noting that the Weyl ordering is 
invariant under the similar transformations, we obtain another generalized differential expression related 
to the Hermite polynomials. Those identites can be applied to investigate the nonclasscial properties of 
quantum optical fields. 

Keywords: Weyl mapping rule; Hermite polynomials; similar transformation 
PACS numbers: 42.50.Dv, 03.67. a, 42.50.Ex 


1 Introduction 


As the ’’language” of quantum mechanics, Dirac’s bra-kets have come to represent a quantum world of 
abstract ideas and universal concepts and can get a far better understanding of quantum mechanics. Because 
conceptions in quantum mechanics quite differ from those in classical mechanics, it is inevitable that quantum 
mechanics must have its own mathematical symbols which are endowed with special physical meaning. For 
instance, in math [1], an inhomogeneous Fredholm equation (FF) of the first kind is written as g (t) = 
f-a ^ (b s) / (s) ds, K (t, s) is the continuous kernel function. In quantum mechanics [H |3], we introduced 


an operator Fredholm equation defined as G (a, = jR ^ (®> of 9) (?) dq, in which the kernel function 

K (a, ab q) is a quantum operator, g is a real variable, a and denote the annihilation and creation operator 
of a quantized radiation field. As is well known, integrations over the operators of type |) (| cannot be directly 
performed by Newton-Leibniz integration rule. 

In Ref. [1], Fan proposed the technique of integration within an ordered product (IWOP) of operators 
which enables Newton-Leibniz integration rules directly working for Dirac’s ket-bra operators with contin¬ 
uum variables. The technique of IWOP shows that the operator Fredholm equation (OFF) can directly per¬ 
form integration if K (a, ob g) is an ordered product operator. An example of taking K (a, ab g) =: exp (g — Q) 
we have 



(.q-Qf 


:/(g)=:G(Q):, 


Q 


a -b 


( 1 ) 


where : exp — (g — Q)^ : is the integral kernel, commutes with a within the normally ordered symbol 

exp — (g — Q)^ : = |g) (g| and the completeness relation ^'^^dq\q) (g| = 1, we 


” : : ”. Noting that 
can see 


v^- 



poo poo 

{q-Qf ■f{q)= dq\q){q\f{q)= dg |g) (g| / (Q) = / (Q), 

J —oo j — oo 


( 2 ) 


where |g) = tt exp (—g^/2-b •\/2ga'^ — a'^^/2) |0) denotes the coordinate representation with its eigen¬ 
function Q |g) = g |g), and : G {Q) : is the normally ordered expansion of / [Q). When / (g) = i7„ (g), we 
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have 



{q-Qf : = 2": 


(3) 


where 2": Q": is the normally ordered form of the operator Hermite polynomials {Q). (q) is the 

single-variable Hermite polynomials with its generation function 


^ n 

X! (^) = (2^9 - 9^) (4) 

n=0 


and Hn (q) spans an orthonormal and complete function space, namely -^e Hn (q) Hm {q) = 2"nlS„m- 
Noting the annihilation operator a = (Q + iP) /v^, and {q\ = 2“^/^ (^q — (q| , we can derive the matrix 

element of {Q) in the coordinate representation and the vacuum state 


(qI (Q) |0) = 2-(g| : : |0), 


and then it follows 


(g) = 2”/2e«'/2 {q\ : (at + a)” : |0) = {q\ a^” |0) = 


= -9V2 


= (-l)"e« ^ 

^ ' dq^ 


(5) 


which is just the differential expression of single-variable Hermite polynomials. The above derivation shows 
that Dirac’s symbol and its own arithmetic rule can promote the development of the basic quantum theory. 
In Ref. [5], author introduced a two-variable Hermite polynomials in complex space 


m.in(m,n) 

Hm,nia,a*)= ^ 


j.lr).! 


1=0 


l\ (to — /)! (n — l)\ 


/ 1 m—l *n—l I 

(—1) a a , a = q ip^ 


whose generating function is 


— r-rHm,n (a, a*) = exp {-tr + ta + ra *). 

?T? rj. 


m,n=0 


( 6 ) 

(7) 


Two-variable Hermite polynomials can be applied in many fields of physics. For instance, Hm,n {o^,ck*) is 
proved to be the eigenmode of the complex fractional Fourier transform laiziii], so it may be observed 
in the light propagation in graded index (GRIN) medium, and this eigenmode is also the mechanism for 
two-dimensional Talbot effect demonstrated in GRIN medium [9]. 

Due to there exists some similarities between the generation function of single-variable Hermite polyno¬ 
mials Hn {x) and that of two-variable complex Hermite polynomials Hm,n (of, a*) , it is interesting to see that 
the differential expression of Hm,n (o:, a*) is similar to that of Hn (x). In this work, by virtue of Weyl’s map¬ 
ping rule and quantum operator Hermite polynomials, we obtain the differential expression of Hm,n (ctyCX*) 
and another generalized forms. So our work is arranged as follows. In Sec. 2, we briefly introduce Weyl 
correspondence rule, which is related to the Weyl ordering and the technique of integration within Weyl or¬ 
dered product (IWWP) of operators. We reveal that quantum correspondence operator of classical function 
/ (p, q) can directly be obtianed by replacing q and p in / (p, q) by Q and P with the function form invariant. 
In Sec. 3 we introduce the two-variable Hermite function. By formulizing the Weyl correspondence, we 
derive a differential expression of Hm,n (ci, o^*). Noting that the Weyl ordering is invariant under the similar 
transformations, the generalized differential expressions are obtained in Sec. 4. Enlightened by those new 
identities, in Sec. V we investigate the nonclassical properties of an excited squeezed vacuum state, which 
can be generated in quantum systems with restricted dimensions. 


2 Weyl Ordering and Weyl Correspondence Rule 

As is well known, the Weyl correspondence rule [TOl [11] , i.e. 

F{P,Q)= j J dqdpf (p, q) A (g, p), (8) 
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is a recipe for quantizing a classical function / (j>, q) defined in classical phase space as a quantum correspond 
operator F {P, Q). A {q,p) is the Wigner operator difined as 


A (p, q) 



(9) 


When / (p, q) = q"^p'', from Eq. (j8|) its quantum corresponding operators is 


q p 


ml 


2J 


qm-lpVQl ^ 


ml 


: Q^-^P^Q‘ : = : 




( 10 ) 


in which : : denotes the Weyl ordering symbol, and Eq. (US tell us that the Weyl ordered correspond 

operator of classcial function / (p, q) is obtained by just replacing p, g in / (p, q) by P, Q with the function 
form invariant. As one of the definite operator orderings (such as normal ordering, anti-normal ordering and 

Weyl ordering), Weyle ordering is a useful one and within Weyl ordering symbol : : the Bose operators 

are permutable. Enlighted by the technique of IWOP, m proposed the technique of integration within Weyl 
ordered product (IWWP) of operators. From (O we obtain the Weyl ordering form of the Wigner opertor 

A(p,g)= :6(j)-P)6{q-Q) :. (11) 


Therefore one can easily obtain quantum correspond operator of classical function h (p, q) by replacing q ^ Q, 
p ^ P, i.e. 

F[P,Q)= : hiP,Q) : = / / dpdqh {p, q) A {p, q). (12) 

J J —OO 

Eq. ([T|) can tell us that the Weyl correspondence rule is also an OFE. Noting that 


Tr[A(pi,gi) A(p 2 ,g 2 )] 


j ~ (^1 A {P2, <?2)] \z) 

^5{qi -g2)<5(pi -P 2 ), 


it then follows that the reciprocal relation of the Weyl correspondence rule is 


(13) 


2TTTr [F(P,Q)A(g,p)] =27rTr 


J J dgidpi/i(pi,gi) A(pi,gi) A(p,g) 


= h (p, q). 


(14) 


In many cases, taking a = {q + ip) /v^, the Wigner operator in can be rewritten as 

d^z 


A(p,g) ^ A(a,a*) = 


\a -\- z) {a — z\ e' 


,a.z —cx z 


= — : exp [—2 (a^ — a*) (a — a)] ' “ 2 ' ^ ~ S {a — a) 

It then follows that the Weyl correspondence formula in Eq. (HU can be recast to 

G(a,a^)= : /(a,a^) : =2 J J dqdpf {a,a*) A {a,a*), 

with its reciprocal relation 


27rTr [G (a, a^) A (a, a*)] = 27rTr : / (a, a^) : A (a, a*) 


= / (a, a*). 


(15) 


(16) 


(17) 


Especially, when G (a, a^) = p (a,a^) describes a density of states for an interesting quantum system, from 
the above reciprocal relation we can see 


27rTr [p (a, a^) A (a, a*)] = W (a, a*), (18) 

W (a, a*) denotes a quasi-probability distribution Wigner function. Thus Eq. (HU can also be called the 
Wigner-Weyl correspondence rule. 
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3 New Differential Expression of Two-variable Hermite Polyno¬ 
mials (a, a*) 

In quantum theory, Hm,n ict,c(*) is the generalized Bargmann representation of the two-mode Fock state in 
the bipartite entangled state representation [T3], i.e. 


|m,n) = ^=r |0,0) 


y/mW. 'Jm\n\ 

and it spans an orthonormal and complete function space, 

■12 




J J (^2^, V2r) = Vmlnlm'ln'lSm,m'Sn,n'■ (19) 


Eg. dT^ indicates any one function / (a, a*) can be expanded by those orthogonal basis, 

OO 

f{a,a*)= (^a, V2a*) . 


( 20 ) 


m,n—0 


where Cm,n is a constant to be determined by follow derivation. From its generation function shown in 0, 
we can also expand the normally ordered form of Wigner operator in Eq. (na as 


1 1 

A (a, a*) = — : exp [—2 (ai — a*) {a — a)] : = —^ 


V2’"+"ai™a” 


),!n! 


Hr 


m,n—0 


(V2a,V2a*^ : . ( 21 ) 


Substituting (OIH) and (1^ into (ITB)) . we have 
^ ^ In ^ m!n! 

7n,n—0 


Hm,n (V2a, V2a*) : ^ C^',n'H:^,{V2a, V2a*) 


m' ,n'=0 


OO OO 


— 2: ^ X/ ^rn',r 

m,n=0 m' ,n'—Q 


V2’"+"ai™a” f (Pa 


i\n\ 


(^a, '/2a*) (V2a, 72^2^) 


= : ^ Cm,u'//^/'^a'^-. =:F{a/a):. 

m,n=0 

Taking the coherent state expectation values of (1^^ . we have 


(a|:F(ai,a): |a) = F (a*, a) = (a| : ^ : |a) = ^ C,^,nV2^ {a*) 

m,n—0 m,n—0 


m n 

a , 


and then 


1 


dm Qn 


(J — _ 

™’" TO!n!\/2™+" da*™- da'' 


F{a*,a) 


Q:=0 


Therefore, from (EH)) we obtain 


/ \ I 

F(a.a-) = j: 


m,n—0 


(23) 


01=0 


This is a new formula for deriving Weyl’s classical correspondence of normally ordered quantum operators. 
For example, when Gi ,a) =: Fi (ai,a) : =: ai^a": , from Eq. (051) we have 


gi (a, a*) = ^ Hm,n ('/2a, V2a*) . 

And then considering the reciprocal relation in (1171) . we can see 

gi {a, a*) = 27rTr [Gi [a/aj A (a, a*)] = 27rTr [ai’”a"A {a, a*)] = „ (V2a, y/^a*) 

V ^ 


(24) 
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On the other hand, given any operator F (a,a^) , its classical correspondence can also be obtained by 
utilizing the following formular (see Appendix for details) 


Substituting Gi (a’l',a) =: : into (1^ . we have 

gi {a, a*) = J ^ (-^| Gi (a^ a) |/3) e2(“/5*-«*/3) = g 2 |a|= J f£ ^ 

= g2|apy' ^ ^_^*^m^ng-2|/3|2+2(a/3*-a*/3)^ 

Making the substitution \/2/3 —>■ /3, —>■ /?*, the above equation can be derived as 

-\P\^ + V2iap* -a*P) 


(25) 


-PI : a1'™a”: IS) 


9iia,a ) = J —P P exp 

(_l)m-»g2|a|= Qm Qn r ^2^ 

- 2™+" 5a™ 5a*"' J tt 

Utilizing the following integral formula 


• exp 


+ V2 {aP* - a*P) 


■ exp 


1 


= -exp 


we can obtain 


91 {a, a*) = 


h |a|^ + sa + 77 a 


2™+" 5a™ 5a* 

Comparing Eq. (l26l) and (l24t . and setting -y^a —>■ a, -y^a* —> a* we can derive 

, f)m f) 

H la a*! - —_- 

which is just the differential form for the generation function of two-variable Hermite polynomials, and for 
the especial case of to = n 


ST]’ 

h . 

, Re [h] < 0, 


- exp 

(-2|a|>)^ 

(26) 

we 

can derive 


,n ( l«l ) ’ 

(27) 


Hr 


r,m {a, a*) = exp (|a|") exp (- |a|") . 


4 Generalized Differential Expressions related to Hermite Poly¬ 
nomials 

In order to generate non-symmetric quantum mechanical representation, in Ref.|18j authors proposed a 
non-unitary operator U, defined as 


U = 


J_ [fl 

y/J^J tt 
1 

exp 

r/JI I 2g 


a T j \ z 


z 

* 


——aP^ + f — 1 ] a^a + 


a 

2fi 


(28) 


which is a quantum operator image of the classical symplectic transformation (z, z*) —>■ {pz + vz* ,az + tz*) 


in phase space, and where 


z 

z 


= exp (za'I’ — z*a) |0) denotes the coherent state [19] and 
^ ^ ^ ~ ~ invers of IJ reads as 


a T I \ z 




z 

* 


T —U 

—a ^ 


2 

z 




: exp 


-a^^ + 


t ^2 

a'a - a 

[2t 


2t j 


(29) 
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(30) 


From (1^51) and (P^ . we can find W ^ and IJ engenders a similar transformation 

UaU ~^ = IJ,a + Ua ^ U ~^ = aa + ra ^ 


and its invers transformation 


U ^aU = Ta — va^, U = ua'^ — 


fia ’ — aa , 


(31) 


where four complex parameters satisfies — va = \ for keeping [^a + va\ aa + ra^] = l.Also it is impor¬ 
tant to note that the Weyl ordering is invariant under the similar transformations [llinillS], i.e. 


lJG{a\a)U ^ = IJ \ g{a\a)j \ U ^ = : <? (era-I-ra^,/xa-I-: . 


(32) 


4.1 Generalized Differential Expression Related to the Product of Two Single¬ 
variable Hermite Polynomials 

Supposing an operator G 2 (a^,a) = |0) (0| and from Eqs. (IT71) and (I5T]) . we can express the 

classical Weyl correspondence of G 2 { a ^, a ) as 


g 2 {a*, a) = 27rTr [G 2 (a^ a) A (a, a*)] = 27rTr a^'^t/ |0) (0| {a, a*) 


= 27rTr 

By virtue of the following formular 


U — aa) 10) (0| (ra — U ^A(a, a*) 


(fa + ga ^) = -i 




we can see 


(/ia'I'- era)"* = ■Hm(^ 

(ra-jea^)" = (^- 


a I g . 
—a — \ —a' 
2g V 2er 


„ a — \ l — a ^ 
2 v V 2t 


Therefore we have 

[—aa + fia^) 10) (0| (ra — 

= {\IW) {]Iy) ^ 

where |0) (0| =: exp [—a^a] : . Utilizing the following formular (for details in Appendix) 

F(^a\a)= : /(a^,a) : =2 / - : (—z| F (a^,a) |z) exp [2 (az* — a^z + a^a)] 

we can obtain the Weyl ordering form of (—era + ^af)"* |0) (0| (ra — rea'^)", namely 
{-aa + /xa^’)™ |0) (0| (ra - zea^')" = 2 (“\/t) 




— Izf — 2a^ z + 2az* + 2a^c 


An explicit integral formula proved in Ref. 


min m,n 


2‘^^m\n\ 


(33) 


(34) 

(35) 


(36) 


(37) 


(38) 


l\ (m — l)\ (n — l)\ 


Hm-i{X*)Hr,-i{X) (39) 
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can be utilized for the derivation of Eq. ([ 55|) . Furtherly, we have 


/ 2 min[m,n| 

—-Hm (az*) Hn il3z) exp [- {z - A) (z* - A*)] = 

^ 1=0 


( 4 q:/ 3 )* 




/! (m — ly. {n — ly. 


Hm-i (aX*) Hn-i il3X), (40) 


then it follows 


/ I -\ / I —\ ’ 

(—(TO + /io^)”^ |0) (0| (to — ^o’^)" = 2TO!n! \\ J ’^ ' 


X : 




U (m — ly. {n — ly. 


Hm-I { V — Hn-i ( \l—a ) exp [-2a'^a] : . (41) 


Noting that Weyl ordering is invariant under the similar transformations shown in Eq. (I32[) . and using the 
Weyl correspondence formula in Eq. (13, we can obtain the classical correspondence of G 2 (o^, o) , 


52 (£t',a) = 2irTr[Gj (ot.o) A(a,Q')l =2irTr [ot"(/|0)(0|(/“'a"A(a,Q") 


= 2m\n\ — 


X Hri.— l 


n min m.n 


E 


Hm—l 


2r ^ 

— [Iia + ly'a j 
V 


l\ (m — iy.[n— ly. 

exp [—2 (era -|- ra*) {fxa + i^a*)]. 


{a a + Ta*) 


(42) 


From Eq. (1^ . the classical correspondence of G 2 {a^, a) can also be expressed as 
92 {a*, a) = ^(-/3|a^™17|0) (0| 17-^0” 1/3) 


= ^e2|“lV^(-/3*)™/3"exp 

y/W J TT 

e 2 |“l" (- 1 )’"+’" r epp 

2‘m+n Qa^ 9a*" J tt 

And then using the following integral formula 


exp 


- 2a*P + 2aP* - —/3" - —P* 


- 2a*p + 2ap* - —P^ - —P* 


(43) 


ePe 


■ exp 


h |a|^ -I- sa -I- 77a* -I- fa^ + ga* 


- ^fg 


exp 


—hsrj -I- s^g + 77^/ 
- 4.fg 


whose convergence conditions are Re{h ± f ± g) <0 and Re ^ ^ < 0, we have 


/ -| \m+n am an 

r * X 2iQ(in-ij o a 

02 (ta , a) = e ' ' —-—^-7;-7;-exp 

^ ^ ^ 2""+” da”^ 9a*" 

Comparing Eqs. (1^ and (lUl) . we can see 

ani C5i 

4 /iT \af + 2 vTa*‘^ + 2 ag,a^' 


—Agr \a\ — 2vTa*^ — 2aga^ 


(44) 


exp 

= 2{2ga)~ {2 tv)^ 


. 9a'" 9a*" 


exp 


—Agr |a|^ — 2vTa*^ — 2aga^ 


min m.n 


Z=0 


X Hm . — ] 


— (era + ra*) 


I \l 




II -y AML 


(45) 


{ga + va*) 


which is a generalized differential formula related to the product of two single-variable Hermite polynomials. 
If taking t = g*, a = ly* to meet \g\'^ — = 1, U is unitary. Thus Eq. (HSl) reduces to 


exp 4 |/7|^ |a|^-I-2/7*72a*^-I-2/77e*a^ 

min[m, 

= 2{2g,y*)^ {2g*iyy ^ 


9a'" 9a*" 


exp 


—4 \g\^ |a|^ — 2g*iya*‘^ — 2gv*OL^ 


;=o 


X Hrn. — l 


2li / * I * * \ 

— (v a-\- ^ a ) 


m\ j n 

i)[i 

Hr,.-] 


n - 


2\g\ 


( 46 ) 


^ (ga + iya*) 


7 















































Especially, for m = n, Eq. is 


exp 


4 1^1^ \a\^ + 2iJ*va*'^ + 2fii/*a 


.*„2 


= 2™+iforfor 


E 

1=0 


m\ Im 

l)\l 


Qm Qm 

. da^ da^'^ 

i 


exp 


-4 \n\^ \a\^ - 2^jL*va*^ - 2nv*a 
1 ^ 


*^2 


in- 


2\fi\ 




Hr,. 


(47) 


4.2 Generalized Differential Expression Related to Two-variable Hermite Poly¬ 
nomials 

Following we shall derive another gereralized differential expression related to two-variable Hermite Polynomi¬ 
als. From Eqs. (IMl) and (I55|) we can express the classical correspondence of G 2 (a^ a) = |0) (0| U~^a^ as 

32 (a*, a) = 27rTr (l — aa) |0) (0| (ra — 17“^A(a, a*) 


= 27t{-J^] I-W^I Tr 


UHra ( J |0) (0| Hr, ( J (a, a*) 


Due to 


we can obtain 


[W2] / .xfc I 

Hm (X) = (2,)-- 




A:! (m — 2k)\ 


32 (a*, a) = 27r(^- 

xTr 




(-t) 



^ to kW{m-2k)lll^{n-2iy. 

U\m-2k){n-2l\U-^A{a,a*) . 


(48) 


Supposed a projection operator of the number state \m) {n\ and from Eq. dSll), its Weyl ordering reads as 

d^z ■ 


\m) {n\ = 2 

= 2 


(—^1 to) (n \z) exp [2 (az* — z + a^a)] 


^2 2; • ^ z'^ 


J TT • 

By virtue of the following integral formula 


• exp 


—2 |z| + 2 (az^ — a' z + a^a) 


(49) 


^^O^knl 

- p p exp 

TT 


-h |/3|' + s/3 + //3*] = (-z)"+' h-^-^e^Hk,i ^ 


we have 


) (n| = : Hm,n (2a^, 2a) exp (—2a^a) : . 


V uIto! 

Therefore, considering Eqs. (l32l) and dsni) , we can see 

32(a ,a) 47r( ^ E E fci _ 2fc)!/! (n - 21)! 


(50) 


k=o 1=0 


xTr 


Hm-2k,n-2i [2 (cTo + Ttt^) , 2 (/itt + i^c^)] exp [—2 (tjo + To^) pa + : A (a, a*) 

Weyl correspondence rule in (fT71) tells us that the classical correspondence of G 2 (a'l', a) is 
32 (a*, a) = 2 (“i) exp[-2(cra-f TQ!*)(/rQ!-t:^a*)] 


ml (-^)'n! 

^ ^ k\{m-2k)\l\in-2l)\ 
A;=0 Z=0 ^ ^ ^ ’ 


Hm-2k,n-2i [2 (aa + ra*), 2 (^a + z/a"')]. (51) 


X 



































Comparing Eqs. m and (HU) , we also derive a simplified equation 


exp 


4/xr \a\^ + 2 i'Ta*^ + 2 afj.a^ 


Qm Qn 

da^ da*^ 


exp 


—4/xt \a^ — 2 vtq*^ — 2cr/ra^ 


[m/2][n/2] m\ ('_JLVn' 

= ^ ^ TTFZT —gm- 2 fc,»- 2 i [2 (era + ra*), 2 (^a + z/a*)] . 


fc=0 i=o 


k\{m-2k)\l\{n-2l)V 


( 52 ) 


The right of (15^ is a summation of two-variable Hermite Polynomials, while Eq. (HSl) is related to the 
product of two single-variable Hermite Polynomials. Eqs. (HSl) and (15^ are new generalized differential 
expressions related to the Hermite Polynomials. For a special case of m = n, Eq. (15^ reduces to 


exp 


4^r |a|^ -I- 2vTa*'^ + 2cr/ia 


i] S’" S’" 

. Sa’" Sa*™ 


exp 


—4/xt |a|^ — 2vTa*'^ — 2afia‘^ 




(-|l) 


k,l—0 


k\ {m-2k)\l\{m-2l)\ 


Hm-2k,m-2i [2 {( 7 a + Ta*) , 2 {pa + ua*)]. 


For a special case of unitary operator, e.g. tj = exp [| — a^)], we have r = /r* = coshr, a = v* = sinhr, 

Eq. (15^ can be simplified to be 


exp 


1 3 ^ 3 ^ 

4cosh^ r lar -I- 2sinhr coshr {a^ -I- a*^) —-—-exp 

' ’ J Sa*" Sa*" 


—4 cosh^ r |a|^ — 2 sinhr coshr (a^ -I- a*^) 


[m/2] [n/2] / tanhr\^+^ 


= 2 cosh' 


m+n ^ 


EE 

A;=0 /=0 


^_ tanh r ^ ^ 


7.lr? J 


k\{m-2k)\ l\{n-2l)\ 


Hm- 2 k,n- 2 i [2 (a* cosh r -|- a sinhr), 2 (a coshr -I- a* sinhr)]. 


5 Applications of Eq. (l45l) in the Study of Nonclassical Features of 
Quantum Light Field 

For a composite system consisting of a multi-level atom (or quantum dot) coupled to a cavity and driven by 
a weak coherent field, quantum-optical effects can be demonstrated in the interaction processes of photon 
emission and absorbtion with atom between ground and excited states [inunjiiiiig. In Ref. [23] , author has 
considered a weak coherent incident field, and the quantum state of emitted light can be expressed as a series 
of excited coherent states Fm 1*^) ^ superposition of the different Fock states j')/') = Fn |n) 

(due to |a) = Fn l’^))- Here we consider a strong coupling case, and the initial state of the incident 
source is a single-mode squeezed vacuum field. Only considering an effective measurement to optical field 
(e.g. photon-statistics methods), this multi-photon processes originated from quantum nonlinearity can be 
monitored by 

p (r, n) = C-^a^-S (r) |0) (0| (r) a". (53) 

which can exhibit similar behavior to that of an excited quantum state (e.g. a^’" ]</?)). Cn = Tr[p{r,n)] = 
n! cosh" rP„ (coshr) is a normalized constant. Pn (coshr) is the expression of the Legendre polynomials. 
S (r) = exp (a^^ “ denotes a unitary operator with the following transformation identities 

S~^ {r) aS {r) = a cosh r-|-sinh r, 

S~^{r)a'S{r) = a’^ coshr-j-asinhr. (54) 

In order to investigate the non-classical features of p (r, n) , Eq. (fT51) shows that its quasi-probability distri¬ 
bution Wigner function can be written as 


IT (a, a*) = 27rTr [p (r, n) A (a, a*)] 


= 271(7" ^Tr 


S' (r) (a^ coshr-I-a sinhr)" jo) (01 (a coshr-I-sinhr)" S' ^(r)A(a,a*) (55) 
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Comparing (IMl) with (EU, we can see ^ = r —> coshr, u = ;/—>■ — sinhr. From Eqs. (1551) and (1331) . thus we 
can obtain the Wigner function of p (r, n), 


W{a,a*) = 


Pn (cosh r) 


— sinh r 


exp 


—2 lasinhr — a* coshrl 


E 

i=0 


n \ 2* (— cothr)* 


I J {n — l)l 




tanhr 


(a* cosh r — a sinh r) 


In Fig. 1, negative part in certain region of phase-space indicates an evidence of nonclassicality of the 
state generated by adding photon to a weak squeezed radiation field. The variance of photon-addition 
can exhibit different nonclassical features with the fixed squeezing value of r = 0.2, which dominates the 
quadrature distribution in the directions of —X and —Y. The top row with three pictures shows that even 
photon have been added to the weak squeezed radiation held, odd photon at bottom row. 
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7 Appendix: Derivation of Eq. (I37l) 


The Glauber-Sudarshan P representation, as one of most important quantum phase space theory, is the 
quasiprobability distribution in which observables are expressed in normal order. By virtue of the coherent 
state representation \z) = exp [za^ — z*a) |0), the P-representation of a quantum density matrix p is dehned 
as 

^= y 1^) (^1 ■ (56) 

The inverse relation of Eq. dnn) is 

P (z, z*) = J yy (-/3| p 1/3) exp (^|/3|^ + z/]* - z* , (57) 

which was hrst obtained by Mehta [25]. Substituting (EH) into (IMl) and Utilizing Weyl ordering of the 
coherent state projector, namely 

|z) (z| = 2 : exp [—2 (z* — a'^) (z — a)] : , 

we have 

p = (-/3|p|/3)exp(^|/3|^-kz/3*-z*/3) : exp [-2 (z* - ot) (z - a)] : 

= 2 y —— : (—/3| p 1/3) exp [2 (a/3* — a^/3-I-a^a)] : . (58) 

which can conveniently recast operators into their Weyl ordering. 
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n = l, r = 0.2 


n=3, r = 0.2 


n = S, r = 0.2 


Figure 1: Wigner distributions of the excited squeezed vacuum state with fixed squeezing. In the top row, even 
number photons have been added to the squeezing va^^um field. The case of odd number photons has been 
shown in the bottom row. 



